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Abstract 

By using the density matrix renormahzation group method, we investigate ground- and excited-state properties 
of the Bg-orbital degenerate Hubbard model at quarter filling for two kinds of lattices, zigzag chain and ladder. In 
the zigzag chain, the system is effectively regarded as a decoupled double chain of the S—1/2 antiferromagnetic 
Heisenberg model, and the spin gap is approximately zero, similar to the case of weakly coupled Heisenberg chains. 
On the other hand, in the ladder, the spin correlation on the rung remains robust and the spin gap exists. 
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Frustration among competing interactions is a key is- 
sue to bring richness of cooperative phenomena in mag- 
netic systems [1]. For instance, antiferromagnets with 
the triangle-base structure possess frustration, since it 
is impossible to satisfy all the interactions. To suppress 
the efi'ect of frustration, spins compromise with each 
other to take the noncoUinear 120° structure. On the 
other hand, when electrons partially fill degenerate or- 
bitals, it is necessary to consider the orbital as an inter- 
nal degree of freedom in addition to the spin. The spa- 
tial anisotropy of the orbital causes electron hopping 
between orbitals depending on the direction, hence it 
is important to take account of the lattice geometry to 
understand the spin-orbital structure even in a simple 
cubic lattice. In particular, it is interesting to consider 
the effect of the orbital anisotropy on the frustrated 
lattice. 

In this paper, we investigate the eg-orbital degen- 
erate Hubbard model at quarter filling on the zigzag 
chain and the ladder. When the Hund's rule coupling 
J is small, it is found that the ground state is a param- 
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agnetic (PM) phase with ferro-orbital (FO) ordering. 
Here we focus on the PM phase and set J—0 for sim- 
plicity. The effect of J will be discussed elsewhere in 
the future. Then, the Hamiltonian is given by 



i,a,7,7' ,a 

+ U ^Pi7TPi7i +U' Piaapiba 



(1) 



where diaa {diba) is the annihilation operator for an 
electron with the spin a in the ci322_^2 {d^2_y2 ) orbital 
at the site i, pi-yrT=d\^^dir.i„, and a is the vector connect- 
ing adjacent sites. Note that the lattices are arranged in 
the (i, i/)-plane as shown in Fig. 1. The hopping ampli- 
tudes are given by tL=i/4, V3t/4, t?6=3t/4 
for the a;-direction, tL=t/4, tL=*L=^*/4, ^66=3^/4 
for the j/-direction, and i"a=i/4, ia6=i"a=\/3t/8, 
f"b=3f/16 along the zigzag path. Hereafter, t is taken 
as the energy unit. We employ the density matrix 
renormalization group method with the finite-system 
algorithm [2], which allows us to investigate the frus- 
trated system with high accuracy. Here we consider 
each orbital as a virtual site to reduce the size of the 
superblock Hilbelt space. The number of states kept 
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Fig. 1. Lattice configurations and numbering of the site. The 
zigzag chain and the ladder are considered as two-chain systems 
with the zigzag path and the rung, respectively. 

for each block m is up to m=200 and the truncation 
error is estimated to be 10~^ at most. 

In the strong-coupling region, each site is occupied 
by one electron, indicating that only spin and orbital 
degrees of freedom remain effective. In Fig. 2(a), the 
spin correlation function Cspin(i, j) = {<S'f ) with Sf — 
'll-,{Pii'\~Piil) /"^ is shown for U=U' =20. In the zigzag 
chain, we find that the antiferromagnetic (AFM) cor- 
relation exists between intra-chain sites in each of a 
double chain, while the spin correlation between inter- 
chain sites is much weak. The double-chain spin struc- 
ture is understood from the FO structure. In order to 
characterize the orbital shape, it is useful to define an 
angle 6i at each site as 



tan0i = {T^-)/{rO 



(2) 



in the spirit of the mean-field approximation [3] , where 
71"=E.(rfLrfi^<.+<.'^i-)/2 and T.^=T.MaAa.- 
d\f^^diba) /i. By using 6i, the new operators are defined 
as Jiaa=e'®'''^[cos(6li/2)diaCT+sin(Si/2)di6a] and di6a = 
e*'''/^[-sin(6li/2)di„<,-fcos(6li/2)di6^] [4], which deter- 
mine the orbital shape. For the zigzag chain, Si~47r/3 



is observed, indicating the (3a; 



-type FO structure 



shown in Fig. 2(b). The orbital shape extends along 
the double chain, not along the zigzag path, to gain the 
kinetic energy by electron hopping. Namely, the AFM 
superexchange interaction along the ti-direction Ji is 
much weaker than that along the a;-direction J2 . Note 



that J2/ Ji is estimated to be 64 for the (3a; 



-type 



FO structure. Thus, due to the orbital anisotropy, the 
spin correlation along the zigzag path is weakened to 
suppress the effect of frustration. On the other hand, 
the spin correlation on the rung remains robust in the 
ladder, since the orbital shape extends to both the leg- 
and rung-directions, as shown in Fig. 2(c). 

In order to consider how the orbital state is affected 
by the spin excitation, we investigate the orbital shape 
in the lowest-energy state with Sfot^l, where Sfot is 
the z-component of total spin. In the zigzag chain, the 
FO structure with the same orbital shape as that in 
the ground state is observed. Namely, the double-chain 
nature remains even in the spin-excited state. Thus, 
it is expected that the spin excitation should be de- 
scribed by the zigzag spin chain with large J2/J1. For 
the zigzag spin chain, it has been found that the spin 
gap decreases exponentially with increasing J2/J1 [5]. 
In Fig. 2(d) , the system-size dependence of the spin gap 
is shown. We find that the spin gap converges to almost 
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Fig. 2. (a) The spin correlation functions for tlie lower- and 
upper-chains from the center of the lower chain in the PM 
ground state. The FO structures for (b) the zigzag chain and 
(c) the ladder, (d) The spin gap vs. 1/iVx with N:c=N /2. The 
dotted line denotes the linear extrapolation to N^ — oo. 

zero in the thermodynamic limit, as expected by anal- 
ogy with the zigzag spin chain. On the other hand, the 
FO structure occurs also in the ladder, although the 
orbital shape slightly changes from that in the ground 
state. A finite spin gap is observed in the ladder, as 
shown in Fig. 2(d). Details of the spin-orbital structure 
in the spin-excited state will be discussed elsewhere. 

In summary, the zigzag chain is effectively decou- 
pled to two chains by orbital ordering to suppress the 
effect of frustration. Accordingly, the spin gap is ap- 
proximately zero in the zigzag chain, while a finite spin 
gap exists in the ladder. 
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